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We analyze the tunneling of vortex states from elliptically shaped traps. Using the hydrodynamic 
representation of the Gross-Pitaevskii (Nonlinear Schrdinger) equation, we derive analytically and 
demonstrate numerically a novel type of quantum fluid flow: a jet-like singularity formed by the 
interaction between the vortex and the nonhomogenous field. For strongly elongated traps, the 
ellipticity overwhelms the circular rotation, resulting in the ejection of field in narrow, well-defined 
directions. These jets can also be understood as a formation of caustics since they correspond to a 
convergence of trajectories starting from the top of the potential barrier and meeting at a certain 
point on the exit line. They will appear in any coherent wave system with angular momentum and 
non-circular symmetry, such as superfiuids, Bose-Einstein condensates, and light. 
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Topological charges, such as vortices, are fundamental 
to the dynamics of coherent fieldsii^. They appear in 
laser systems, carry charge in superconductors, charac- 
terize turbulence in quantum fluids, and hold potential 
for quantum memory"^. To date, the main focus in vor- 
tex dynamics has been on transport, so that the charges 
could move and interact, (see e.g. Ref. |4|) However, it 
is often desirable, and sometimes necessary, to confine 
and trap vortex structures. This is a basic problem in 
trapping theory, yet it has received very little attention. 
Here, we consider the dynamics of vortex decay in a po- 
tential and show that asymmetry in the potential can 
lead to the development of jets during wave tunneling. 
These formations concentrate wave density in the form 
of caustics and represent a new type of coherent structure 
for wave transport. 

The emphasis on vorticity implies that phase dynam- 
ics will be important to the tunneling process. Even 
in the context of simple wavefunctions, without angu- 
lar momentum, phase can have profound effects. Ex- 
amples include the recent prediction of "blips" in the 
outgoing matter through a trap^""~ and the develop- 
ment of dispersive shock waves^rJ^, e.g. when tunneling 
through a barrier ;i2, These latter structures are traveling 
waves with oscillating phase that are finding increasing 
importance in fluids^ii^, optics^iii, and Bose-Einstein 
condensatesi^iii. In spatially inhomogeneous potentials, 
such as the elliptical wells typical of BEC experiments^^, 
both shock waves and blips can go unstable and generate 
vortices. Here, we consider the simplest case of a circu- 
lar vortex trapped in an elliptical well and examine the 
competition of symmetry during wavefunction tunneling. 

Tunneling problems are usually discussed within the 
framework of the WKB approximation, which looks 
for a solution of the Schrodinger equation (Nonlinear 




FIG. 1: A cross-section of the potential Vtrap{x, y) of an ellip- 
tical trap with cylindrically symmetric internal potential well, 
which holds a vortex state shown symbolically. The potential 
Vweii (x, y), which tends to a constant value for large x and/or 
y (shown by dashed-doted lines) is used for the preparation 
of the initial state. Vb is the difference between the top rim of 
the barrier and the chemical potential level /x whose crossing 
with the barrier determines the exit line. The insert shows a 
view from above with the circular top rim, a vortex inside it, 
and the elliptic exit line outside. 



Schrodinger (NLS) or Gross-Pitaevskii (GP) equation in 
our case) 

ih-^ = - — y^^l^ + Vtrap{v)^ + \m^^P (1) 

as an expansion over h. In the case of the NLS equa- 
tion describing propagation of a classical coherent elec- 
tromagnetic wave h — 1, ip is the amplitude of the elec- 
tric component, the propagation distance z plays part of 
time, the part of mass m is played by the wave vector 
of the light emitted by the laser, and the potential V 
is created by changing the linear refraction index of the 
medium. The appearance of singularities near the turn- 



ing points is typical for WKB approach. These singu- 
larities make generalization to 2d and higher dimensions 
a formidable technical problem (see, e.g. Ref. [l3 and 
references therein) since then the points become singular 
lines or surfaces. (Nevertheless calculations can be quite 
straightforward in the cylindrical symmetry case^^ even 
from the vortex state.) The most promising way to deal 
with this problem is to use the hydrodynamic approach, 
which was applied to the tunneling dynamics in Refs. 
I5li2ll [See also Refs. I9l-[l2|. Two hydrodynamic equations 
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follow directly from the GP equation ([T]) for the density 
p = l^p and the velocity mv — —hVip fields determined 
by the phase (p of the wave function ^. The principal 
advantage of this approach is the absence of any singular 
behavior of the wave function near the turning points. 
For example, the temporal tunneling dynamics was de- 
duced analytically^i^ from these equations accounting 
for the role of the interaction (nonlinearity) in Id sys- 
tems. 

Here we consider wave tunneling from a two- 
dimensional trap. We first consider an irrotational initial 
state in order to present our approach and for the sake 
of comparison. We then generalize it to the rotational 
case with a vortex initial state. We will show that in the 
case of nonzero angular momentum, jets appear in which 
matter/intensity exit the trap along certain preferential 
directions. 

In the Id case^i (see also a more detailed derivation in 
Ref. 5), the adiabatic approximation in the Euler equa- 
tion ([2]) yields the velocity at the exit point in the form 
v{xex) — v'2Vo/to, where Vq is the energy difference be- 
tween the top of the barrier and the chemical potential 
in the trap. The same definition holds for the 2d case 
(see Fig. [1}. In this higher dimension, there is an exit 
curve where the chemical potential p crosses the trap po- 
tential. Integrating the exit flux /9(r)ve2:(r) • fi over this 
closed curve (n is the unit vector normal to this curve), 
gives the total exit flux 



dN{t) 
dt 



- f p{r){-Vexir) ■n)dl. 



(4) 



The 2d tunneling problem^ is mapped onto the classical 
motion of a fluid droplet (tracer) falling down from the 
top of the difference potential AV{x,y) = Vtrapix^y) — 
Vweiiix^y) to the observation points {x,y) on the exit 
line. Here, Vweii{x,y) is the potential well used to pre- 
pare the initial state and Vtrap{x, y) is the actual poten- 
tial of the trap from which the tunneling takes place (see 
Fig. [T|) . In the case of an irrotational flow, the exit veloc- 
ity vector is found from the equation mv = — VAy(r). 



Considering the example of the elliptic trap shown in Fig. 



^V{x,y) 



-(2/^+eV) 
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with the aspect ratio e < 1 we get w^x{x,y) — [ecux^coy). 
Adiabatically slow varying density^i^i in the trap may 
be described as p{N,r,t) = N{t)po exp{—^^/2mVo\r\) 
and Eq. ^ becomes 



dN 
'dt 



= ~NL, 



(6) 



where lirrot is the integral escape rate. Using the polar 
coordinates, (r, 0), the differential escape rate reads 



dlirrot 2Vb sin2 61 -F £3 cos2 61 --4f= 
dO muj x^i^T^) 



(7) 



where Au = Vo/{huj) and x(£;^) = sin ^ -|- e^ cos^ 0. 
This equation is obtained by using the simple connection 
eta.n/3 = tan6' between the polar angle 9 and the angle 
/3 of the exit velocity direction (see Vex after Eq. 5). 

Assuming certain dependence (say, linear) of the in- 
teraction on the number of particles N remaining within 
the trap, its dependence on time for various values of 
the parameters of the system can be found from Eq. 
dl]). Qualitatively these dependencies are rather close to 
those considered in detail in Ref. ^] , discussing tunnel- 
ing from one dimensional traps, although the numerical 
values may differ. 

A 2d configuration allows the consideration of a ro- 
tational initial state, i.e. one with nonzero vorticity 
(D = V X V, which is not possible in Id. Since the ve- 
locity field is TTiv — —K^Lp, a finite vorticity in Eq.Q 
appears if the phase Lp is singular along some lines in 
3d space or points in 2d space. As the wave function 
is single valued, the equation m f v(r)dl = 2T:hiy with 
an integer v holds for integration over any closed path. 
The phase ip = — ^^^ arctan - corresponds to a vor- 
tex around the line x = y = with the velocity field 
v™t(r) = iy:^{~y,x,0) where g = ^yx^ +y'^. The vor- 
ticity reads l3 — i'—hz5{x)5{y). Eq. (|3]) requires that 
the density field vanishes us p{q) cx g^^"^ at p — >■ and 
the quantum potential blows up as g""^ . 

The total velocity field v(r) — Vrot(r)+Vpot(r) is a sum 
of the rotational and irrotational velocity fields with V x 
"^pot — 0. Then we may carry out the same program as 
above, i.e. we prepare the initial state in a potential well 
yweiiix,y) with a vorticity characterized by an integer 
quantum number v. Then we change the potential into 
Vtrap{x,y) allowing the wave function to tunnel through 
the barrier and apply the adiabatic approximation.i^^ 
As a result, the Euler equation ([3]) takes the form 

■m — '—^^ + yV[vpof(r,i) -I- v^ot(r)]^ = 



-V 



AF(r) - 



2m (p- 



(8) 



with the centrifugal potential in the right hand side. Now 
the Cole - Hopf transformation!, allows one to map the 
tunneling problem on the classical motion of a fluid tracer 
described by the equation 



mwrot + rnwpot — wv x a; — V 



Ay(r) - 



2mg'^ 



(9) 



The "Lorentz force" in (|9]) is zero everywhere except for 
the line x = y = and does not play a role. 

For -Vpot — and Ay(r) — 0, one gets -Vrot = ^ • 
fiz X ri-i{t) with Trotit) = g {cos fit, sin nt,0) and il — 
— ^, which corresponds to the tracer making a circular 
rotation with the velocity v = hv/mg. We again use 
the elliptic potential ([5]) and consider the simplest vortex 
with V — 1] then the escape rate is calculated similarly 
to the rotationless case. We are interested in the angular 
dependence of the differential escape rate 



dL. 



dp 




e cos /3 cos 6* + sin /3 sin 9 



de 

dp 



Vx(i 



(10) 



Contrary to the irrotational case (O, there is not now a 
simple relation between the polar coordinate 9 and the 
angle /3 of the escape direction. Therefore, we have to 
use the more general equation (jlOp in which the depen- 
dence I3{9) is found by solving equation of motion (|9]) 
numerically. 

The angular dependence ([TUl) of the escape rate from 
a vortex state is the principal result of this paper. It 
shows that the rate is determined largely by the func- 
tion /3(0), which is not necessarily monotonous. Fig. [5] 
shows that at small enough aspect ratio e and moder- 
ate barrier height Vq, the escape direction (3(9) obtains a 
maximum value Pmax on the exit curve at an angle 9max- 
At this point, the derivative jg blows up, indicating that 
an ensemble of streamlines originating from a spread of 
9 angles around 9max collapses together. This bunching 
results in the formation of a jet in the /3max direction. 

For small eccentricity and moderate vorticity v, the 
exit flow is distorted only weakly. For large enough ec- 
centricity and low barrier heights, or high enough vortic- 
ity, this distortion becomes strong and jets appear. An 
example of well-developed jets are shown in Fig. |31 ob- 
tained when the shape of the trap strongly deviates from 
the circular one (e = 0.1). 

The jets are caused by the interplay between the el- 
liptical shape of the trap and the spherical symmetry of 
the vortex. If the equivalent tracer motion were domi- 
nated by the circular motion only, then we could have 
a strange situation in which the tracer could have left 
the elongated trap, traveled along its circular orbit, and 
then tried to re-enter the trap. Interestingly, this sce- 
nario is prevented by the irrotational part of the velocity 




FIG. 2: Exit angle /3 as a function of polar angle 6 for the bar- 
rier hight Am varying from 0.5 to 4 for a highly eccentric trap 
(aspect ratio e = 0.1). All the curves pass maxima indicating 
the jets coming out of the trap. 
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FIG. 3: Normalized polar graphs of the differential escape 
rate ^ as a function of /3 for e = 0.1. 



field, whose contribution leads to a frustration point and 
the formation of the jets. The meeting of many differ- 
ent tracer trajectories in the vicinity of 9max on the exit 
line gives rise to the formation of caustics. The resulting 
interference is similar to the tunneling dynamics consid- 
ered in Ref. [19|. As shown in Fig. |31 the variation in 
tunneling speeds gives rise to wave steepening and a new 
type of angular shock. 

In principle, the exponentially weak tunneling from an 
elliptically shaped trap violates the spherical symmetry 
of the problem and may cause a decay of the vortex state. 
In practice, however, the time span of the decay can 
be rather long (esp. if the vortex can be stabilized by 
a strong enough interaction)^^, and observation of the 
above dynamics should be possible. Currently, the best 
candidate systems to observe jet-like tunneling are cold 
atoms in elliptically trapped Bose-Einstein condensates 
and coherent light confined in elongated optical waveg- 
uides. 





FIG. 4: Numerical calculation of the evolution of a light propagating in an elliptically shaped semitransparent waveguide for 
the vorticity u = 1 (left) for An/no = 5.0 • 10"^ written in PMMA and i/ = 8 for An/no = 1.0 ■ 10~*(right) induced in SBN. 
Formation of jets with the propagation of light at a distance 7.5 mm is clearly seen. 



Here, we numerically demonstrate angular caustic for- 
mation in the optical case. The cylindrical waveguide 
required for an experimental realization can be written 
inside a medium such as a glass or polymethyl methacry- 
late (PMMA) using femtosecond laser pulses. As a real- 
istic example, we consider a waveguide with a transverse 
refractive index profile (potential well barrier) shown in 
Fig. 1. The inner side of the well is circular while the 
outer edge is an ellipse with a sem i-major axis of AO/jLir 
(corresponding to an eccentricity Vl — e^ = 0.9). Simu- 
lations were carried out for the 2-|-ld system by solving 
the NLS equation ([l} using a split-step beam propagation 
code. Results are shown in Fig. 21 Even with a small 
refractive index difference (potential barrier height) of 
An/uQ = 5.0 • 10~^, jet-like tunneling occurs for a singly- 
charged vortex v = 1 (see Fig. 21 left panel). For compar- 
ison, we also give similar results for light propagating in 
an optically-induced waveguide in SBN (Strontium Bar- 
ium Niobate). In this case, the writing beam diffracts. 



so the waveguide diameter is limited. Even with a 50/im 
inner radius, the potential still diffracts. This weaker po- 
tential means that we cannot see the jet form if the input 
vortex is only singly-charged. Shown is the output field 
of a charge ^ = 8 vortex after propagating 8mm in the 
crystal (see Fig. JH right panel). 

In summary, by using the hydrodynamic formulation of 
the nonlinear Schrodinger equation, we were able to carry 
out an analysis of tunneling from a vortex field trapped in 
an asymmetric potential well. Interference between the 
rotational motion of the field with a strongly asymmet- 
ric tunneling flow created angular caustics, resulting in 
jet-like radiation patterns. Analytic results were verified 
with numerical simulation. 
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